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In this paper, from the microscopic Hubbard Hamiltonian we extract the local spin-singlet amplitude 
due to short-range spin correlations, and quantify its strength near half-filling. As a first application 
of the present approach, we study a problem of the energy dispersion and its d-wave modulation in 
the insulating cuprates, Sr2Cu02Cl2 and Ca2Cu02Cl2- Without any adjustable parameters, most 
puzzling issues are naturally and quantitatively explained within the present approach. 

PACS numbers: 71.10.Fd, 71.27.+a 
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Recent discovery of a normal state pseudogap in un- 
derdoped copper oxides jjj has attracted considerable at- 
tention both from experimentalists and theoretical physi- 
cists for many years. For these materials, the low fre- 
quency spectral weight begins to be strongly suppressed 
below some characteristic temperature T* higher than 
T c . This anomalous behavior has been observed through 
various experimental probes such as angle resolved pho- 
toemission spectroscopy (ARPES), |||| specific heat, ||] 
tunneling, 0] NMR, || and optical conductivity. One 
of the most puzzling questions in pseudogap phenomena 
is why T* has a completely different doping dependence 
from T c , in spite of possibly their close relation. 

There are many theoretical scenarios to attempt to un- 
derstand the pseudogap phenomena. These include the 
spinon pair formation without the Bose-Einstein conden- 
sation of holons in the slave boson (or fermion) mean- field 
theory, |@|| strong superconducting (SC) phase fluctua- 
tions, pi strong magnetic fluctuations near the antifer- 
romagnetic (AF) instability, Jll[] and so on. At present 
there is no consensus on the origin of the pseudogap and 
its relationship with the SC long-range order. Among 
several scenarios, the slave boson (or fermion) mean- field 
theory ||[)| may shed some insight into the problem. This 
is because the pseudogap is closely related to a spin gap, 
the predicted phase diagram is, at least, qualitatively 
consistent with experiments, and furthermore it starts 
from the microscopic model (t — J Hamiltonian) as op- 
posed to other phenomenological models. On the other 
hand, the recombination of a spinon and a holon into a 
physical electron is nontrivial, and also the constraint of 
no-double-occupancy at each site is difficult to impose 
at a microscopic level. In deriving a local spin-singlet 
amplitude which may be responsible for the pseudogap 
behavior, in this paper we use the Hubbard Hamiltonian 
instead of the t — J Hamiltonian. Local spin-singlet am- 
plitude is induced directly from short-range spin corre- 
lations in the normal state and its strength is quantified 



near half-filling. In this paper short-range spin correla- 
tion means that when site i is occupied by an electron 
with up-spin (or down-spin), then the nearest sites pre- 
dominantly by electrons with the opposite spin. 

We start by defining the one-band Hubbard model 
proposed by Anderson ||] as the simplest model which 
might capture the correct low energy physics of copper 
oxides. The Hubbard model is described by the Hamil- 
tonian where Ci. a destroys an electron at site i with spin 
a on a two-dimensional square lattice 



H = -t 



E < 

<i,j>,cr 



U ^2 4,^4,1^,1 
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t is a hopping parameter between nearest neighbors < 
i,j > and U denotes local Coulomb repulsion. It is be- 
lieved that the realistic strength of the Coulomb repulsion 
lies in between the weak and strong coupling regimes, 
namely, U ~ W — 2W where W is the bandwidth of 8t 
in two dimensions. 

As a first step to the microscopic understanding of the 
pseudogap, it is important to answer a more fundamen- 
tal question: How can spin-singlet tendency |l2j appear 
directly from the Coulomb repulsion, presumably with- 
out the exchange of bosonic degrees of freedom such as 
spin fluctuations |l5(] ? As a simple example to illustrate 
this point qualitatively, let us consider the U 3> t limit at 
half-filling where only an up-spin or down-spin electron 
is allowed at each site. The typical spin-singlet structure 
(with d-wave form factor) in which we are interested in 
this paper is A g (i) = \ J^s 9{5){c l+&A c lA - c.+^c^), 
where g(S) is an appropriate structure factor. For a d- 
wave symmetry, for instance, 



1/2 if S =(±1,0), 
g(S) = { -1/2 ifa = (0,±l), 
if otherwise. 
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Then the local spin-singlet amplitude for strongly cor- 
related electrons {(\A g {i)\) = (1 X 1 + X 0)/2 = 1/2) 



1 



is increased over its noninteracting value ((|A ff (i)|) = 
(1/2 x 1/2 + 1/2 x l/2)/2 = 1/4 for both spins). 

To make this argument more quantitative, we intro- 
duce a spin-singlet correlation function p4| 



Xs ( i ,r) = (T r A ff (z,r)At(0,0)}, 



(3) 



where T T is the imaginary time ordering operator. In 
this paper we consider only the local spin-singlet am- 
plitude (|A g (0)| 2 ), which may be obtained by Xg(i ~ * 
0,t — > 0~). Now we examine whether there is an in- 
crease in (|A g (0)| 2 ) for strongly correlated electrons with 
respect to (|A ff (0)| 2 )o for the noninteracting electrons, 
which is similar in spirit to a renormalization group (RG) 
approach [jl5| . As already illustrated in the previous 
paragraph, the increase in (|A ff (0)| 2 ) for strongly cor- 
related electrons crucially depends on the short-range 
spin correlations. Although there exists no controlled 
way of obtaining the wave function in general, certain lo- 
cal or short-range static quantities such as double occu- 
pancy (or equivalently local spin amplitude) or the near- 
est neighbor correlations are reasonably well captured by 
the mean-field state with AF order. Before going further, 
it is important to establish the validity of the current ap- 
proximation by explicitly comparing the above quantities 
calculated in the mean-field state (with AF order) with 
available quantum Monte Carlo (QMC) results. 

The double occupancy D = (n^n^i) plays an impor- 
tant role in gauging the degree of strong correlations in 
the Hubbard-type model. For instance, D — > at half- 
filling for U — > oo, while Dq evaluated in the noninteract- 
ing state is (n/2) 2 = 0.25. In Fig. |l](a) our calculations 
(solid curve) are compared with (virtually exact) QMC 
results (open circles) by White et al. [|l6| for n = 1, 
JV = 4 x 4, and T — t/16. For purely interaction in- 
duced effect, Do is subtracted from D. Above U ~ it 
the agreement with QMC data is increasingly better and 
for U > 8t the two results are almost indistinguishable. 
For the nearest neighbor correlations between i and j, 
we calculate t e ff/t = {c\ a Cj^)u I (4,<r c J>V=o for com- 
parison with available QMC results. Again t e ft/t eval- 
uated in the noninteracting state, namely, unity is sub- 
tracted from t e ff/t for purely interaction induced effect. 
Figure 0(b) shows our results evaluated in the mean- 
field state (solid curve) together with QMC results (open 
circles) by White et al. |L6| From intermediate to very 
strong coupling, the deviation with QMC results is less 
than 10% level. In fact the interaction induced local spin- 
singlet amplitude (|A g (0)| 2 ) -(|A s (0)| 2 )o involves mainly 
nearest neighbor static correlations between electrons. 
Note that the spin density wave (SDW) approximation 
will be used below only to capture the reasonable local 
and short-range correlations between electrons. 

We should also point out some limitations in this ap- 
proach. The physical reason that the nearest neighbor 
correlation is overestimated (by up to 10%) compared 



with QMC data is as follows. Suppose site i is occu- 
pied by an electron with up-spin in the strong coupling 
limit. Then the nearest sites are occupied by down-spin 
electrons as majority and also by up-spin electrons as mi- 
nority. As the distance from site i increases, the effective 
polarization of electron spins decreases in magnitude. In 
the mean-field state with AF order, however, the effec- 
tive polarization is the same for all the distances from 
site i. This means that physical quantities evaluated in 
the mean-field state become progressively overestimated 
for increasingly distant correlations. Another limitation 
is that the mean-field critical doping x c where AF mean- 
field order vanishes increases with increasing U. This is 
inconsistent with more advanced treatment of the Hub- 
bard model [fj"7f where with increasing U, x c decreases 
after reaching its maximum around at x ~ 0.21. Hence 
our approximation based on the Hubbard model is valid 
near half-filling (probably up to a; ~ 0.1 — 0.15). For 
the t — J Hamiltonian the current approximation will be 
valid even far away from half-filling, because this strong 
coupling feature is already taken into account there. 

Then for d-wave type symmetries, the local spin-singlet 
amplitude (|A ff (0)| 2 ) evaluated in the mean-state with 
AF order becomes 

" £ #(*)[/(25- (fc)) + f(E + (k))] - sign[^ {l ~ Qh 



4N 



2UN^" K ' A(fc) 



£^(fc)-f^[/(£-(fc)) - f(E+(k))} , (4) 



where 



A(fc) = yJ((e$)-e(U+Q))/2) 
E±(k) = (e(fc)+e(fc + Q))/2±A(fc) 



2 _1_ A 2 

sdw ' 
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s(k) = — 2i(cos k x + cos k y )— /z for nearest neighbor hop- 
ping, [i is the chemical potential controlling the particle 
density n, N the total number of lattice sites, f(E) the 
Fermi-Dirac distribution function, Q the AF wave vector 
(tt, 7r), 4> g {k) the Fourier transform of gs(i), and the sum- 
mation accompanied by the prime symbol is over wave 
vectors in half of the first Brillouin zone. The SDW gap 
A s( i w and the chemical potential fi are self-consistently 
determined through the gap and number equations for 
given U, T and n. Figure||(a) shows the local spin-singlet 
amplitude (|Ad| 2 ) subtracted by that evaluated in the 
noninteracting state for purely interaction induced effect. 
The most interesting finding is that the local spin-singlet 
amplitude is induced directly by short-range spin corre- 
lations without any explicit driving (attractive) interac- 
tion. The short-range spin correlations are ultimately as- 
cribed to the strong local Coulomb repulsion U. The in- 
teraction induced (|Ad| 2 ) rapidly decreases with increas- 
ing doping, and if the second limitation in the current 
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approximation mentioned above is properly corrected, it 
will decrease even faster with doping. 

Because the value of the first term is nearly the same as 
(|A+(0)| 2 )o for the noninteracting electrons, the positiv- 
ity of the second term [p| leads to an increase in the local 
spin-singlet amplitude. For d xy {<j){k) = 2 sin k x sin k y ) 
symmetry, for example, the local spin-singlet ampli- 
tude is even suppressed due to <j){k + Q) = <j){k). We 
found that among several symmetries including local 
s-wave, extended s-wave, d xy , and d x 2_ y 2, the d x i_ y i 
(<Mfc) = cosfc x cosfcy) symmetry shows the largest 
increase in the local spin-singlet amplitude. Likewise the 
local spin amplitude ((S z (i)) 2 ) is also induced from the 
same Coulomb repulsion. 

In order to determine the effective strength of the in- 
duced local spin-singlet amplitude (with d-wave symme- 
try) corresponding to Fig. ||(a), we consider a Hamilto- 
nian 



H = -t J2 </M + ^E A dW Ad W 



(6) 



with the same parameters (i, T, and n) as before. Now 
the strategy is to find Vmd in such a manner that Vi n d 
gives the same interaction induced local spin-singlet am- 
plitude (|A(i| 2 ) -(|A,j| 2 }o as in Fig. |2](a) . Please note that 
Vmd is the effective strength of the induced local spin- 
singlet amplitude, but not attractive interaction yield- 
ing a local SC pair in the Bose-Einstein limit. As be- 
fore, we may gauge the validity of a given approximation 
(BCS approximation in this case) by explicitly comparing 
some local and short-range static quantities calculated in 
that approximation with available numerical data for the 
above Hamiltonian. Unfortunately there exist no such 
numerical results. Then one may use general knowledge 
which has been obtained through study on the repulsive 
and attractive Hubbard models in the SDW and BCS 
approximations, respectively. The SDW approximation 
for the repulsive Hubbard model which was discussed in 
details in the previous paragraphs, suggests that for inter- 
mediate to strong coupling the mean-field approximation 
reasonably well captures both local and short-range static 
quantities. For the BCS approximation of the attractive 
Hubbard model 19 1, it has also been shown that the local 
static quantity (double occupancy) in that approxima- 
tion is in good agreement with QMC data for \U\ > 3t. 
However, the numerical results of nearest neighbor cor- 
relations in the attractive Hubbard model are not avail- 
able in the literature to our knowledge. In spite of lack 
of complete numerical data for the attractive Hubbard 
model and for the Hamiltonian Eq. |^, one may safely 
expect that the BCS mean-field approximation to Eq. ^ 
will give at least reasonable local and short-range static 
quantities for intermediate to strong coupling. Note that 
the mean-field (BCS) approximation will be used below 
only to capture the reasonable local and short-range cor- 



relations between electrons, as before in the SDW ap- 
proximation. 

In the BCS approximation for the above Hamiltonian, 
(|A d (0)| 2 ) is given as 



n 

AN 



^f) tanh (£W/2T)] 



[^Mk)^MuME{k)/2T)Y , (7) 



where E{k) = \J e 2 {k) + A 2 ,(fe). As before, the d-wave 

gap A d (k) — A(f> d (k) and the chemical potential ju are 
self-consistently determined through the gap and num- 
ber equations for given Vi n d, T and n. In Fig. 0(b) Vi n d 
is plotted as a function of doping x = 1 — n for T = 
and U = 1.5W = I2t near half-filling. The BCS ap- 
proximation (| Vind | > 3£) is expected to be accurate for 
x < 0.1 and to be less accurate beyond it. The strength 
of the induced local spin-singlet amplitude rapidly de- 
creases with doping, just as the induced AF correlations 
do. Beyond x ~ 0.1 — 0.15, \Vi n d\ will decrease even faster 
when the strong coupling effect is correctly included, as 
was mentioned in Fig. ||(a). For (7 > t at n = 1, Vi n d 
saturates to be — 7.64i and (|A,j(0)| 2 } — ► 0.25, consistent 
with our previous result ((|Ad(0)|) = 0.5) based on qual- 
itative argument. Vi n d and its doping dependence have 
some interesting consequences which may answer some 
of the puzzling issues in the high temperature supercon- 
ductors. 

As a first application of the present approach, we 
attack a problem of the energy dispersion and its d- 
wave modulation in the insulating cuprates, Sr2Cu02Cl2 
and Ca2Cu02Cl2- Recent ARPES experiments for an 
insulating cuprate S^CuC^C^ ^(j clearly show that 
the near isotropy and the overall band dispersion along 
(7r/2,7r/2) - (tt,0) and (7r/2,7r/2) - (0,0) cannot be ex- 
plained by considering only AF order or its fluctuations, 
unless some adjustable fitting parameters such as t' and 
t" are introduced. Furthermore a d-wave-like modulation 
of the insulating gap in Ca2Cu02Ci2 |2l[] is totally myste- 
rious from that point of view. Because the induced local 
spin-singlet and spin amplitudes increase with decreasing 
doping, near half-filling and at low temperatures, both 
short-range spin-singlet and AF fluctuations are strong 
and coexist, leading to a resonating valence bond (RVB) 
type state M. Our general experience |E| tells that as 
long as the energy dispersion is concerned, it is basically 
identical in a strongly fluctuating state (or in the pseu- 
dogap state) and in a long-range ordered state. The only 
difference is that in the former state the gap is filled with 
some spectral weight (pseudogap) and the spectral func- 
tion is spilt into two (relatively broad) peaks instead of 
two delta functions. Thus at half-filling and at low tem- 
peratures, a hole strongly interacts with both short-range 
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spin-singlet and AF fluctuations, yielding an energy dis- 
persion similar to that in the coexistence state of the SC 
and AF long-range order. p|| The energy dispersion is 
then given by 

^e 2 (k) + A 2 dw + A 2 (cosk x ~cosk y ) 2 . (8) 

From the overall bandwidth (~ 320 meV) experimen- 
tally measured from (tt/2,tt/2) to (0,0) along which 
4>d(k) vanishes, the Coulomb repulsion U (and A sc i w — 
U/2) is determined to be 11 At for t =250 meV. For 
U = 11 At, n = 1 and T = 0, the effective strength 
Vi n d of the induced spin-singlet amplitude is completely 
determined to give Vind = — 6.09£ and A = 2.07i. 
Figure || clearly demonstrates that the energy disper- 
sion and its near isotropy along (ir/2,ir/2) — (tt,0) and 
(7r/2,7r/2) — (0,0) are in excellent agreement with ex- 
periments without any adjustable parameters. Our en- 
ergy dispersion along (n/2, ir/2) — (n, 0) is proportional to 
(cos k x — cos k y ) 2 (solid curve in the inset) at low energies 
as opposed to | cos k x — cos ky | (dashed curve) predicted 
from the flux-phase in the t — J model. p3] Recent ex- 
periments show a significant deviation of the energy dis- 
persion from the cusp-like form along (tt/2, ir/2) — (n, 0). 
A similar result to ours was recently obtained in the con- 
text of SO (5) symmetry. |25|] The state with a spin gap 
(or pseudogap) of (A(j)d) 2 /U ~ Jcf) 2 d and also with a Mott- 
Hubbard gap of order U at half-filling is expected to con- 
tinuously evolve into a state with a relatively smaller spin 
gap (or pseudogap) but without a charge gap away from 
half-filling. 

In this paper we have considered only the local spin- 
singlet amplitude induced by short-range spin correla- 
tions and its consequences. The long-range ci-wave su- 
perconductivity is beyond the scope of the present ap- 
proach. How local spin-singlets acquire local SC phases 
and eventually establish their long-range phase coherence 
is a challenging problem to the theory of high tempera- 
ture superconductivity. Quantitative analytical study of 
the Hubbard model in the physically relevant regime and 
of the interplay between AF and ci-wave pairing correla- 
tions beyond a mean- field level is not available at present, 
mainly due to the absence of a small parameter. Cer- 
tainly it is a subject for future study. It is also desirable 
to perform numerical simulations for a Hamiltonian with 
ci-wave symmetry such as Eq. |s| if" it is possible. Then the 
mean-field BCS approximation used here to obtain rea- 
sonable local and short-range static correlations between 
electrons will be tested for its range of validity. 

In summary, from the microscopic Hubbard Hamilto- 
nian, the local spin-singlet amplitude due to short-range 
spin correlations was explicitly extracted and its strength 
was quantified near half-filling. As a first application of 
the present approach, a problem of the energy dispersion 
and its d-wave modulation in the insulating cuprates, 
Sr 2 Cu02Ci2 and Ca 2 Cu02Cl2 was studied. Without 



any adjustable parameters, most puzzling issues on the 
energy dispersion were naturally and quantitatively ex- 
plained within the present approach. 
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FIG. 1. (a) Double occupancy D — (rii^ni^i) and (b) near- 
est neighbor correlation t e f f /t = (c- itr Cj jCT )t//(4 iCT .Cj, CT )t7=o cal- 
culated in the mean-field state with AF order (solid curve) and 
in the QMC simulations by White et al. jl6| (open circles) for 
n = 1, N = 4 X 4, and T = t/16. For purely interaction in- 
duced effect, Do and unity are subtracted from D and t e ff/t, 
respectively. 




FIG. 2. (a) Local spin-singlet amplitude (lA^I 2 ) calcu- 
lated in the mean-field state with AF order for U = 12t, 
T = 0. (|Atj| 2 )o is subtracted for purely interaction induced 
effect, (b) Effective strength Vind of the local spin-singlet am- 
plitude induced by short-range spin correlations for U = 12t, 
T = 0. 
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FIG. 3. Calculated energy dispersion for an insulating 
cuprate Ca2Cu02Ci2. The dispersion is measured from 
(tt/2, 7r/2) point with A sd „ = U/2 and A = 2.07t. The solid 
and dashed curves in the inset denote the dispersions calcu- 
lated from our approach, and from the flux phase in the t — J 
model, respectively. 
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